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Abstract 

We present measurements of the azimuthal dependence of charged jet production in central and semi¬ 
central = 2.76 TeV Pb-Pb collisions with respect to the second harmonic event plane, quanti¬ 
fied as vf' Jet finding is performed employing the anti-A:T algorithm with a resolution parameter 
R = 0.2 using charged tracks from the ALICE tracking system. The contribution of the azimuthal 
anisotropy of the underlying event is taken into account event-by-event. The remaining (statisti¬ 
cal) region-to-region fluctuations are removed on an ensemble basis by unfolding the jet spectra for 
different event plane orientations independently. Significant non-zero is observed in semi¬ 
central collisions (30-50% centrality) for 20 < p'^ < 90 GeV/c. The azimuthal dependence of the 

charged jet production is similar to the dependence observed for jets comprising both charged and 
neutral fragments, and compatible with measurements of the V 2 of single charged particles at high 
Pt- Good agreement between the data and predictions from JEWEL, an event generator simulating 
parton shower evolution in the presence of a dense QCD medium, is found in semi-central collisions. 
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1 Introduction 


The aim of the heavy-ion program at the LHC is to study strongly interacting matter in ultra-relativistic 
nuclear collisions where the formation of a quark-gluon plasma (QGP), a deconfined state of quarks and 
gluons, is expected 12]. Hard partons that propagate through the collision medium lose energy via (mul¬ 
tiple) scattering and gluon radiation lii. Jet measurements are used to experimentally explore parton 
energy loss in the hot and dense medium. Studies at the LHC and RHIC have shown that jet and high-px 
single particle production in heavy-ion collisions are suppressed with respect to the expected production 
in a superposition of independent pp collisions lIMl. This observation is consistent with energy loss, 
which is further supported by measurements of dijet energy asymmetry and di-hadron angular correla¬ 
tions lEM]. 

In non-central Pb-Pb collisions, the initial overlap region of the colliding nuclei projected into the plane 
perpendicular to the beam direction has an approximately elliptic shape. Jets emitted along the minor 
axis of the ellipse (defined as the in-plane direction) on average traverse less medium - and are therefore 
expected to lose less energy - than jets that are emitted along the major axis of the ellipse (the out- 
of-plane direction). The dependence of jet production on the angle relative to the second-harmonic 
symmetry plane ^2 (the symmetry plane angles *P„ define the orientations of the symmetry axes of the 
initial nucleon distribution of the collision) can be used to probe the path-length dependence of jet energy 
loss. This dependence is quantified by the parameter v'^ the coefficient of the second term in a Louder 
expansion of the azimuthal distribution of jets relative to symmetry planes 'P„, 


dN 
h ( (Ojet 


oc 1 -f ^ 2d„'^‘cos [n ((pjet - %)] , 

n—l 


( 1 ) 


where (pjet denotes the azimuthal angle of the jet. 

In central collisions, the average distance that a jet propagates through the medium is approximately 
equal in the in-plane and out-of-plane directions, therefore a small is expected. In semi-central 
collisions the average in-medium distance is shorter, while the relative difference between the average 
distances in-plane and out-of-plane is larger, hence a non-zero is expected. Lluctuations in the 
initial distribution of nucleons within the overlap region can lead to additional contributions to and 
higher harmonic coefficients in the Louder decomposition. 


The path-length dependence of parton energy loss is of particular interest because it is sensitive to the 
underlying energy-loss mechanism. Lor collisional (elastic) energy loss, the amount of lost energy de¬ 
pends linearly on path leng th, while for radiative (inelastic) energy loss, the dependence is quadratic due 


to interference effects Ill7l .ll8l] . Some strong-interaction models based on the AdS/CLT correspondence 
suggest an even stronger path-length dependence j 19, Earlier studies of the V 2 of high-^^x single 
particles have already tested the path-length dependence of energy loss 121142511. Comparisons of these 
results to theoretical calculations have shown that the V 2 is sensitive to several aspects of the medium 
evolution, including the effects of longitudinal and transverse expansion and the life time of the system 
until freeze-out ilw . It is therefore important to measure multiple observables that are sensitive to the 
path-length dependence of energy loss, such as recoil yields of charged particles and jets jllL 1271 12811. 
Jets are expected to better represent the original parton kinematics and provide more detailed informa¬ 
tion on energy loss. Theoretical predictions from JEWEL, which couples parton shower evolution to the 


presence of a QCD medium with a density derived from Glauber simulations Il29l . 13011 . have shown that a 
finite is expected for non-central collisions at the LHC. Similar results have been found in 142*^' studies 
in heavy-ion collisions generated by the AMPT model |31, 32]. A first measurement of of jets 

comprising both charged and neutral fragments has been reported by the ATLAS collaboration ISSll . The 
results presented in this paper extend the measurement to a lower pj range (px > 30 GeV/c for 
central collisions and pj > 20 GeV/c for semi-central collisions). 
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In this article, measurements of of R = 0.2 charged jets reconstructed with the anti-^j jet finder 
algorithm in Pb-Pb collisions with 0-5% and 30-50% collision centrality are presented. The largest 
experimental challenge in jet analyses in heavy-ion collisions is the separation of the jet signal from the 
background of mostly low-pj particles from the underlying event and from unrelated scatterings that take 
place in the collision. The jet energy is corrected on a jet-by-jet basis using an estimate of the background 
transverse momentum density which takes into account the dominant flow harmonics V 2 and V 3 of the 
background event-by-event, as will be described in Sections IZTl and lZ^ The coefficient is obtained 
from pt- differential jet yields measured with respect to the experimentally accessible event plane *Pep, 2 , 
which is reconstructed at forward rapidities (2.8 < rj <5.1 and —3.7 < rj < —1.7, Sec. l2.1l) . The reported 
V 2 *’ has been corrected back to the azimuthal anisotropy with respect to the underlying symmetry plane 
^^2 by applying an event plane resolution correction (Sec. 12.41) . Jets are reconstructed at mid-rapidity 
(Ifljetl < 0.7) using charged constituent tracks with momenta 0.15 < pj < 100 GeV/c, and are required 
to contain a charged hadron with pj > 3 GeV/c. The in-plane and out-of-plane jet spectra are unfolded 
independently to take into account detector effects and remaining azimuthally-dependent fluctuations 
in the underlying event transverse momentum density (Sec. 12.31) . The jet spectra are corrected back to 
particle-level jets consisting of only primary charged particles from the collision. 


2 Experimental setup and data analysis 


ALICE is a dedicated heavy-ion experiment at the LHC at CERN. A full overview of the detector layout 
and performance can be found in 11331.13411. The central barrel detector system, covering full azimuth, is 
positioned in a solenoidal magnet with a field strength of 0.5 T. It comprises the Inner Tracking System 
(ITS) built from six layers of silicon detectors (the Silicon Pixel, Drift, and Strip Detectors: SPD, SDD 
and SSD) and a Time Projection Chamber (TPC). The two inner layers of the ITS, which comprise the 
SPD, are located at 3.9 and 7.2 cm radial distance from the beam axis. 


The data presented in this paper were recorded in the Pb-Pb data taking periods in 2010 and 2011 at 
= 2.76 TeV, using a minimium-bias trigger (2010) or an online centrality trigger for hadronic 
interactions (2011), which requires a minimum multiplicity in both the VOA and VOC detectors (discs 
of segmented scintillators covering full azimuth and 2.8 < rj < 5.1 and —3.7 < rj < —1.7, respectively). 
The VO detectors are used to determine event centrality based on the energy deposition in the scintillator 
tiles ll35n and the event plane orientation, see Sec. 12.11 Centrality, determined from the sum of the VO 
amplitudes, is expressed as percentiles of the total hadronic cross section, with 0-5% referring to the 
most central (largest multiplicity) events i35n . The trigger is fully efficient in azimuth in the presented 
centrality ranges. Centrality estimation using the VO system does not bias the 'Pep, « determination [? ]. 
Time information from the VO detectors is used to reject beam-gas interactions from the event sample 
and the remaining contribution of such interactions is negligible. Only events with a primary vertex 
position within ±10 cm along the beam direction from the nominal interaction point were used in the 
analysis. A total of 6.8x10^ events with 0-5% centrality and 8.6x10^ events with 30-50% centrality, 
corresponding to integrated luminosities of 18 and 5.6 respectively, are used in this analysis. 


Charged particle tracks in this analysis are measured by the ITS and TPC and are selected in a pseudora¬ 
pidity range |t]| < 0.9 with transverse momenta 0.15 < pj < 100 GeV/c. To ensure a good momentum 
resolution, tracks were required to have at least three hits per track in the ITS. Since the SPD acceptance 
is non-uniform in azimuth for the data sample used in this analysis, two classes of tracks are used. The 
first class requires at least three hits per track in the ITS, with at least one hit per track in the SPD. The 
second class contains tracks without hits in the SPD, in which case the primary interaction vertex is used 
as an additional constraint for the momentum determination. Eor each track, the expected number of 
TPC space points is calculated based on its trajectory; tracks are accepted if they have at least 80% of 
the expected TPC space-points, with a minimum of 70 TPC points. Tracks produced from interactions 
between particles and the detector, as well as tracks originating from weak decays (‘secondary tracks’) 
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are rejected. The contribution of secondary tracks to the track sample is less than 10% for tracks with pT 
< 1 GeV/c and negligible for tracks with higher transverse momentum. 


2.1 Event plane determination 

The coefficient quantifies azimuthal anisotropy with respect to *1^2 ■ The azimuthal anisotropy of 
the underlying event (‘background flow’) is also described by a Fourier series with harmonics v„ = 
(cos(?i[(p —'T„])) 1^ 37] where (p denotes the track azimuthal angle. However, since the initial distribu¬ 
tion of nucleons is not accessible experimentally, the event plane angles 'Tep, the axes of symmetry 
of the density of outgoing particles in the transverse plane, are used in place of when measuring 
and Vn- 

The event plane angles 'Tep, 2 and 'Tep, 3 in this study, corresponding to the two dominant Fourier har¬ 
monics, are reconstructed using the VO detectors. Each VO array consists of four rings in the radial 
direction, with each ring comprising eight cells with the same azimuthal size. The calibrated amplitude 
of the signal in each cell, proportional to the multiplicity incident on the cell, is used as a weight Wceii in 
the construction of the flow vectors Qn OSH 


Qn — VPcellCXp (/ n <Pcell) • 


( 2 ) 

cells 

In order to account for a non-uniform detector response which can generate a bias in the *Tep „ azimuthal 
distribution, the components of the 2„-vectors are adjusted using are-centering procedure |3^ ^]. The 
VGA and VOC detectors cover different rj regions in which multiplicity N and background flow v„ may 
differ. The total VO Q-vector is therefore constructed using weights Xn i38n that are approximately 
proportional to the event plane resolution in each detector. 


2 2 

2n,V0 = Xn,V0AQn,V0A-i- Xn,V0cQn,y0C, 


( 3 ) 


to achieve the optimal combined event plane resolution. The event planes are reconstructed from the real 
and imaginary parts of Q„ as 


1'EP,,=arc.a„(|M)/„, 


( 4 ) 


The V 2 *’ itself is measured with respect to the second harmonic event plane angle. It is corrected for the 
finite precision with which the true symmetry plane is measured in the VO system by applying an event 
plane resolution correction, see Sec. 12.41 


2.2 Jet reconstruction in the presence of background flow 

Jet finding is performed using the FastJet iESi implementation of the infrared and collinear safe kj 
and anti-kj sequential recombination algorithms using the pj recombination scheme and taking massless 
jet constituents. The resolution parameter R = 0.2 determines the characteristic maximum distance of 
constituent tracks to the jet axis in the rj-tp plane. 


In heavy-ion collisions, a large combinatorial background is present from particles that are not related to 
the hard scattering that produced a given jet. This background is subtracted from each jet on an event- 
by-event basis. The anti-kx algorithm is used to find signal jets. A fiducial cut of |T]jet| < 0.7 is applied 
on the signal jets to ensure that all jets are fully contained within the ITS and TPC acceptances and 
edge effects are avoided. The contribution of combinatorial (or ‘fake’) jets (clustered underlying event 
energy) to the measured jet spectrum is reduced by requiring that reconstructed jets contain at least one 
charged particle with pi >3 GeV/c and have an area of at least 0.56 nR^. These selection criteria leave 
the hard part of the jet spectrum unaltered while significantly reducing the number of combinatorial jets 
which stabilizes the unfolding procedure 


m. 
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Fig. 1 : Transverse momentum density of charged tracks as a function of azimuthal angle for a single event from 
the most central 0-5% event class. Data points (blue) are given with statistical uncertainties only. The red curve 
is the fit of Eq. |5]to the distribution, the green and gray curves, obtained from the fit of Eq. |5]as well, show the 
independent contributions of V 2 and V 3 to Peh(^). The dashed magenta line is the normalization constant Pq. 


The /cx-algorithm is used to estimate the average transverse momentum density of the underlying event, 
(Pch)> on an event-by-event basis. The quantity (Pch) is the median of the distribution of P™ dijetM (iho 
ratio of transverse momentum to jet area) of reconstructed R = 0.2 kj-jets, excluding the leading two 
jets from the sample as proposed in |^] and implemented in earlier ALICE jet studies The 

kj jets are required to lie within |T7jet| < 0.7 and have an area A > 0.01. The jet area A is determined 
by embedding a fixed number of near zero-momentum ghost particles per event prior to jet finding; the 
number of ghost particles in each reconstructed jet then gives a direct measure of the jet area. A ghost 
density of 200 particles per unit area is used, so that approximately 25 ghost particles are clustered into 
a jet with a radius of 0.2. 


In each event, the anisotropy of the underlying event is modeled using the dominant ll45ll flow harmonics 
V2 and V3, 

Pch('P) = Po (1 + 2{v2cos [2 (cp - 'Eep. 2 )] + V3 cos [3 ((P - 'Erp, 3)]}) • (5) 


Here, Pch(<P) is the azimuthal distribution of summed track pj for tracks with 0. 15 < px <5 GeV/c 
and Ifltrackl < 0.9. The parameters po and v„ are determined event-by-event from a fit of the right side 
of Eq. [5]to the data. The event plane angles 'Eep, n are not fitted, but fixed to the VO event plane angles. 
A single event example of this procedure is illustrated in Eig. [TJ where the data points represent the 
transverse momentum density distribution in a single event, the red curve represents the full functional 
description of Pch(<P) (Eq.|5ll, the green and gray curves give the contributions of the separate harmonics 
V2 and V3, and the dashed magenta line is the normalization constant po. To reduce the bias of hard jets in 
the estimates of v„ in Eq. [5] while retaining azimuthal uniformity, the leading jet in each event is removed 
by rejecting all tracks for which |T]jet — fitracki < The rj separation between the tracks and the VO 
detectors also removes short range correlations between the event planes and tracks. 


The number of bins to which Eq. [5] is fitted is set on an event-by-event basis to the square root of 
the number of tracks. The fit maximizes the estimated likelihood ll46ll . which is based on a Poisson 
distribution for the bin content. Since the bin contents are not pure counts, but weighted by pj, the 
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Fig. 2: The 5pi distribution (Eq. |8]l from the random cone (RC) procedure as function of cone azimuthal angle 
(f>RC relative to the event plane. In panel (a) the azimuthally-averaged background (pch) has been subtracted; in 
panel (b) the azimuthally dependent Pch(^) from an event-by-event fit of the pj-density with Eq. |5] The solid 
black line represents the mean of the 5pj distribution. 


statistical uncertainties on each bin a, are estimated as the sum of the squares of the px of the individual 
particles: a,- = cj(^pt) = ^ scaled Poisson distribution WiP{xi\}JLi/wi) is used as the probability 

distribution for the data points in the likelihood calculation, with a scale factor w, = of /ji where y,- is 
the bin content and Pi/w, is the expected signal from the fit function. The compatibility of each fit with 
the data is tested by calculating the evaluating the probability of finding a test statistic at least as 

large as the observed one in the distribution. When this probability is less than 0.01, the average event 
background density (Pch) is used instead of Pch(<P); this occurs in 3% (most central) to 7% (semi-central) 
of events. The acceptance criterion is varied in the systematic studies; the sensitivity to it is small. 

The corrected transverse momentum of a jet of area A is calculated from the measured raw jet 
momentum, PjJhjet’ 

pfj'=‘=p-,^,-PchiocalA (6) 

where Pch local is obtained from integration of Pch(<P) around <pjet ±7? 

Pch local — I Pch{(p)d(p. (7) 

Z«Po .J(p-R 

The pre-factor of the integral, is chosen such that integration over the full azimuth yields the av¬ 
erage transverse momentum density (Pch)- The validity of Eq. [5] as a description of the contribution 
of background flow to the underlying event energy is tested by placing cones of radius R = 0.2 at ran¬ 
dom positions (excluding the location of the leading jet) in the rj-cj) plane and subtracting the expected 
summed transverse momentum in a cone from the measured transverse momentum in the cone, 

5pT = Ipr"^-p7r/?'. (8) 

Here, p is the expected transverse momentum density. This procedure is repeated multiple times per 
event, until the full phase space is covered, to obtain a distribution of 5pT values. The Spj distribution as 
a function of the cone azimuthal angle (^rc relative to the event plane 'Tep, 2 is shown in Pig. [2l In panel 
(a) (Pch) has been used for the estimation of the underlying event summed px and in panel (b) Pch(<p). 
Incorporating azimuthal dependence into the underlying event description leads to a sizable reduction in 
the cosine modulation of the 5px distribution. 

The effectiveness of the subtraction of background flow is quantified by comparing the expected and 
measured widths of the dpj distribution in the absence of background flow, a(5px''"^*^), (see Pig. Wp) to 
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Fig. 3: Centrality dependence of the measured and expected relative change in the 5pj disttibution width from 
using the azimuthally dependent Pd, local instead of the median (Pch)- The blue points give the expected reduction 
from simple assumptions about the behavior of charged particle spectra and flow harmonics v„ (following Eq.j^and 
cnii. The red points use the measured widths from 5pj distributions directly. Statistical uncertainties are smaller 
than the marker size. 


the expected and measured widths of the 5pT distribution in the presence of background flow, a(5pT'’") 
(Fig. 12^). Assuming independent particle emission and Poissonian statistics, the expected width of the 
5pT distribution in the absence of background flow (v„ = 0) is given by 1431] 


o{5pl^ °) — \Jnao'^{p'i) +Na{p'iY 


(9) 


where Na is the average expected number of tracks within a cone, (px) the mean pj of a single particle 
spectrum and cj(pt) the standard deviation of this spectrum. This expectation can be extended to include 
contributions from background flow by introducing non-Poissonian density fluctuations (the background 
flow harmonics v„) ll43l] . as 


a(5p]^") = ^Nao'^{pt:) + (Aa + 2A2(v^ + v2))(pt)2. 


( 10 ) 


The measured widths are obtained from the Spx distributions directly; the distributions are constructed 
using as the transverse momentum density p in Eq. [8] either (Pch) to obtain a(5py ) or Pch local for 

o{8pT\ 

Figure [3] shows the expected and measured relative change in the width of the 5px distribution, quan¬ 
tified as (a(5p^') — o{5p'f^'~^))/o{5pj), as function of collision centrality. The blue points give the 
expected reduction from Eqs|9]and[T0l The red points use the measured widths from 5px distributions. 
The expected change is in good quantitative agreement with the measured change over the entire cen¬ 
trality range, indicating that the width of the 5px distributions can be understood in terms of a simple 
independent particle emission model with background flow contributions. 

The background subtraction, unfolding, correction for the reaction plane resolution as described in Sec¬ 
tions 12.31 and 12.41 were also validated using events consisting of PYTHIA jets embedded in heavy-ion 
background events and toy model events. In the first study, full PYTHIA pp events were combined with 
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reconstructed Pb-Pb collisions to create events with a conttolled signal and background. The signal 
jets from PYTHIA have no preferred orientation, = 0, while the heavy-ion events have a non-zero 
V 2 of the soft particles. Jets found in the events were matched to the embedded PYTHIA jets and the 
analysis was carried out with matched jets only. After unfolding, the was compatible with 0, as 
expected. The other study was based on events generated using a simple thermal model for soft particle 
production and a distribution of high-jiT particles that resembles the jet spectrum, as suggested in 14711 . 
A non-zero V 2 = 0.07 was introduced for momenta pj < 5 GeV/c to model the background flow and 
two variations at large pT > 30 GeV/c: V 2 = 0 or V 2 = 0.05. In both cases, the input flow values were 
correctly reconstructed by the analysis. 

2.3 Unfolding 

After the subtraction procedure presented in the previous section, the measured jet spectrum is unfolded 
ll48u49tl to correct for detector effects and fluctuations in the underlying event transverse momentum den¬ 
sity. Mathematically, the unfolded jet spectrum can be derived from the measured spectrum by solving 




( 11 ) 


.gen 


, the unfolded true jet spectrum, where is the measured jet spectrum, 

is a functional description {response function) of distortions due to background fluctuations and detector 
response, and 1^® 1^1 Hiding efficiency. The coefficient is not affected by the effi¬ 

ciency, hence £(p|'^chjet) 1*® omitted from here on. Since the measured jet spectrum is binned, Eq.fTH 
is discretized by replacing the integral by a matrix multiplication 


Mm = G„tTi 


( 12 ) 


where Tj. is the solution of the discretized equation (the prime indicates that Tj. is not corrected for 
jet-finding efficiency). The combined response matrix is the product of the response matrices 

from detector effects and transverse momentum density fluctuations, the latter of which are constructed 
independently for the in-plane and out-of-plane spectra by embedding random cones at specific relative 
azimuth with respect to the event plane (see the text below Eq.[T3]for the definition of the intervals). 


The detector response matrix is obtained by matching pp jets generated by PYTHIA 115011 (‘particle-leveT 
jets) to the same jets after transport through the detector (‘detector-level’ jets) by GEANT3 115IH . where 
the detector conditions are tuned to those of the Pb-Pb data-taking periods. Particle-level jets contain 
only primary charged particles produced by the event generator, which comprise all prompt charged par¬ 
ticles produced in the collision, as well as products of strong and electromagnetic decays, while products 
of weak decays of strange hadrons are rejected. Matching is based on the shortest distance in the rj-cp 
plane between detector level and particle level jets and is bijective, meaning that there is a one-to-one 
correspondence between detector and particle level jets. The response matrix for background fluctuations 
is constructed from the 5pi distributions, which, when normalized, are probability distributions for the 
change of the jet energy caused by background fluctuations. 

Solving Eq. [12] requires inversion of Gjjj j and generally leads to non-physical results which oscillate 
wildly due to the statistical fluctuations of the measured jet yield. The unfolded solution therefore needs 
to be regularized. In general this is done by introducing a penalty term for large local curvatures asso¬ 
ciated with oscillations. Various algorithms for regularized unfolding exist; the unfolding method based 
on the Singular Value Decomposition (SVD unfolding) 15211 is used in this study. A comparison to the 
unfolded solution from minimization 15311 is used in the systematic studies. 


The measured jet spectrum is taken as input for the unfolding routine in the range 30 < < 105 

GeV/c for 0-5% collision centrality and 15 < < 90 GeV/c for 30-50% collision centrality. The 

lower bound corresponds to five times the width of the 5pT distribution, the upper bound is the edge of 
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the last measured bin which contains at least 10 counts. This configuration was found to lead to reliable 
unfolded solutions in Monte Carlo studies Jd, T7]. The unfolded jet spectrum starts at 0 GeV/c to allow 
for feed-in of true jets with low In addition, combinatorial jets which are not rejected by the jet area 
and leading charged particle requirements are migrated to momenta lower than the minimum measured 
^chjet unfolded solution ranges up to 200 GeV/c (0-5%) and 170 GeV/c (30-50%) to allow for 
migration of jets to a higher than the maximum measured momentum. As the data points of the 
unfolded solution are strongly correlated for p^^^^ outside the experimentally measured interval, 
will be reported only within the limits of the measured jet spectra. 


2.4 Evaluation of 

The coefficient is calculated from the difference between the unfolded /jt- differential jet yields in¬ 
plane (Ain) and out-of-plane (Aout) with respect to the second harmonic event plane, corrected for event 
plane resolution, 

chjet, chjet, 71 1 

^2 (Pt ) =- 


^2Ain(p^‘'j'=‘)+A^out(pr^) 

Eq.[T3]is derived by integrating Eq.[T]for n = 2, over intervals [—f, f ] and [^, for Ain and [|, 
and [ 


chjet\ 


(13) 


57C l7t ~\ 
4 ’ 4 


Ml 
4 J 


for Aout^ substituting 'Tep, 2 for *1^2. Eq. [T3]is sensitive to correlations between even-order 
harmonics V2n and 'Tep. 2- As a result of the integration limits however, the first harmonic of the Eourier 


expansion that can contribute to the observed 


IS V. 


. The VO event plane resolution ^^2 is intro¬ 


duced to account for the finite precision with which the true symmetry plane *1^2 is measured in the VO 
system and is defined as 

<%2 = {cO^[2{^ll2-^2)])- (14) 

Measuring event planes in multiple T] regions {sub-events) allows for the evaluation of the resolution 
directly from data i54i 15511. Using the full VO detector and negative and positive r\ sides of the TPC as 
sub-events, the resolution in Eq.[T3]is evaluated as 


cos 


9 I vuvo _ xt/Tpc, I) > 0 
^ ' ^EP, 2 ^EP, 2 


COS 


9 I mvo _ mTPC. ij < o 
^ ' ^EP, 2 ^EP. 2 


52 = 


COS 


9 f mTPC,i,>0 _mTPC.r,<o31 \ 

^l^^EP, 2 ^EP. 2 



1/2 


(15) 


The event plane resolution Si 2 is found to be 0.47 in 0-5% centrality and 0.75 in 30-50% centrality with 
negligible uncertainties. The 'Pep, 2 angles in the TPC are obtained following the procedure of Eq.|4]on 
tracks with 0.15 < px < 4 GeV/c, using unit track weights in the construction of the flow vectors Q 2 
(see Eq. [2]). 

Using the VO detectors for the reconstruction of the event plane guarantees that the jet axis and event 
plane information are separated in pseudorapidity by |At] | > 1 and thus removes autocorrelation biases 
between the signal jets and event plane orientation. The possible non-flow correlation between the event 
plane angle and jets due to di-jets with one jet at mid-rapidity and one jet in the VO acceptance was studied 
using the PYTHIA event generator. The rate of such di-jet configurations was found to be negligible (less 
than 1 per mille of the total di-jet rate at mid-rapidity) for > 20 GeV. Possible effects from back- 
to-back jet pairs with a jet in each of the VO detectors are even smaller. 


2.5 Systematic uncertainties 

The measured is corrected for experimental effects, such as the finite event plane resolution and 
detector effects on the jet energy scale as well as the effects of the uncorrelated background and its fluc¬ 
tuations using the corrections outlined in the Sections 12.1112.41 Hydrodynamic flow of the background 
is taken into account event-by-event in the underlying event description, residual effects are removed 
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by azimuthally independent unfolding. The remaining uncertainties in these correction procedures are 
treated as systematic uncertainties. Systematic uncertainties on v'^ are grouped into two categories, 
shape and correlated, based on their point-to-point correlation. Shape uncertainties are anti-correlated 
between parts of the unfolded spectrum: when the yield in part of the spectrum increases, it decreases 
elsewhere and vice versa. Correlated uncertainties are correlated point-to-point. Both types of uncertain¬ 
ties however have contributions which lead to correlated changes of Ain and Aout- 


Correlated uncertainties are estimated for the in-plane and out-of-plane jet spectra independently. Two 
sources of correlated uncertainties are considered: tracking efficiency and the inclusion of combinatorial 
jets in the measured jet spectrum. The dominant correlated uncertainty (< 10%) arises from tracking and 
is estimated by constructing a detector response matrix with a tracking efficiei^ reduced by 4% (moti¬ 
vated by studies |01 comparing reconstructed tracks to simulations of HIJING |56] events). The observed 
difference between the nominal and modified unfolded solution is faken as a symmefric uncertainly lo 
allow for an over- and underesfimalion of Ihe fracking efficiency. The sensilivily of Ihe unfolded resull 
lo combinatorial jels is lesled by changing Ihe lower range of Ihe unfolded solufion from 0 fo 5 GeV/c, 
which leads fo an overall (correlated) increase of Ihe unfolded jel yield. Bofh correlated uncerfainlies are 
added in quadralure and propagaled to v'^ assuming lhaf varialions are sfrongly correlaled befween 
Ihe in-plane and oul-of-plane jel speclra, while still allowing for effecls from azimulhally-dependenl 
variations in Irack occupancy and reconsfrucfion efficiency, by selling Ihe sample correlation coefficienl 
p = Gij/{OiOj) fo 0.75. 

Shape uncerfainlies fall info Ihree cafegories: assumpfions in Ihe unfolding procedure, feed-in of com¬ 
binatorial jels, and Ihe sensilivily of Ihe unfolded solufion to the shape of the underlying event energy 
distribution. The dominant contribution to the unfolding uncertainty is related to the regularization of 
the unfolded solution. The SVD algorithm 0] regularizes the unfolding by omitting components of the 
measured spectrum for which the singular value is small and which amplify statistical noise in the re¬ 
sult. To explore the sensitivity of the result to the regularization strength, the effective rank of the matrix 
equation that is solved is varied by changing an integer regularization parameter k by ± 1. The SVD un¬ 
folding algorithm uses a prior spectrum as the starting point of the unfolding; the result of the unfolding 
is the ratio between the full spectrum and this prior. The unfolded solution from the algorithm 1531] is 
used as prior (default) as well as a PYTHIA spectrum. The bias from the choice of unfolding algorithm 
itself is tested by comparing the results of the SVD unfolding and the algorithm. 


The same nominal unfolding approach is used for the in-plane and out-of-plane jet spectra and the 
5pT distributions for the in-plane and out-of-plane background fluctuations are similar in width; the 
unfolding uncertainty is therefore strongly correlated between the in-plane and out-of-plane jet spectra. 
These correlations are taken into account by applying the variations in the unfolding procedure to the in¬ 
plane and out-of-plane jet spectra at the same time and calculating the resulting variations of The 

total uncertainty from unfolding is determined by constructing a distribution of all unfolded solutions in 
each interval and assigning the width of this distribution as a systematic uncertainty. 

The other two components of the shape uncertainty are the sensitivity of the unfolded solution to com¬ 
binatorial jets and uncertainties arising from the description of the underlying event; both are estimated 
on the in-plane and out-of-plane jet spectra independently and propagated to as uncorrelated. A 
systematic uncertainty is only assigned when the observed variation is found to be statistically incompat¬ 
ible with the nominal measurement. The effect of combinatorial jets is tested by varying the minimum 
of the measured jet spectrum by ± 5 GeV/c, effectively increasing or decreasing the possible con¬ 
tribution of combinatorial jet yield at low jet momentum. To test the assumptions made in the fitting of 
Eq. |5]the maximum p^i of accepted tracks is lowered to 4 GeV/c. Additionally, the minimum p-value 
that is used as a goodness of fit criterion is changed from 0.01 (the nominal value) to 0.1. The minimum 
required distance of tracks to the leading jet axis in pseudorapidity is enlarged to 0.3. 
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Uncertainty on vf' 


(GeV/c) 

30-40 

60-70 

80-90 

30-40 

60-70 

80-90 


Centrality (%) 


0-5 



30-50 



Unfolding 

0.017 

0.012 

0.016 

0.016 

0.011 

0.015 

Shape 

chiet 1 

Pj -measured 

0.013 

<C stat 

< stat 

0.024 

<C stat 

<C stat 


Pch('P)fit 

0.015 

< stat 

0.016 

<C stat 

<C stat 

<C stat 

Total 


0.027 

0.012 

0.023 

0.029 

0.011 

0.015 

Correlated 

Tracking 

-unfolded 

0.009 

< stat 

0.009 

<C stat 

0.009 

<C stat 

0.007 

< stat 

0.007 

< stat 

0.007 

<C stat 

Total 


0.009 

0.009 

0.009 

0.007 

0.007 

0.007 


Table 1: Systematic uncertainties on Vj*' for various transverse momenta and centralities. Uncertainties in central 
and semi-central collisions are given in the same pj ranges. The definitions of shape uncertainty and correlated 
uncertainty are explained in Sec. 12.51 Fields with the value ‘<C stat’ indicate that no systematic effect can be 
resolved within the statistical limits of the analysis. 

Table [U gives an overview of the systematic uncertainties in terms of absolute uncertainties on for 
all sources (where the total uncertainty is the quadratic sum of the separate components). High statistics 
Monte Carlo testing has been used to verify that uncertainties labeled ‘<C stat’ are indeed negligible 
compared to other uncertainties. 


3 Results and discussion 


The coefficients as function of for 0-5% and 30-50% collision centrality are presented in 
Fig. m Significant positive is observed in semi-central collisions and no (significant) pj depen¬ 
dence is visible. The observed behavior is indicative of path-length-dependent in-medium parton energy 
loss. The observed v'^ in central collisions is of similar magnitude. The systematic uncertainties on the 
measurement however are larger than those on the semi-central data, in particular at lower 
as a result of the larger relative background contribution to the measured jet energy. 


The significance of the results is assessed by calculating a /j-value for the hypothesis that = 0 over 
the presented momentum range. The /j-value is evaluated starting from a modified calculation that 
takes into account both statistical and (correlated) systematic uncertainties, as suggested in 15711 . The 
modified x"^ for the hypothesis = ft, is calculated by minimizing 


X (f^corr) f^shape ) 



iv2,i + ^corr^corr,! “1“ ^shape Mi) 

O? 



e 


2 

shape 


''shape,! 


(16) 


with respect to the systematic shifts £shape> £coit> where V 2 j represent the measured data {n points), a, 
are statistical uncertainties and ashape,i^ O'con-.i denote the two specific types of systematic uncertainties. 
The parameter Cshape is a measure of the fully correlated shifts; a shift of all data points by the correlated 
incertainty acon.i gives a total contribution to x^ of one unit. The systematic shifts for the shape uncer¬ 
tainty are taken to be of equal size for each point, since this gives the best agreement with the = 0 
hypothesis and thus provides a conservative estimate of the significance; the penalty factor is constructed 
such that an average shift of all data points by aghape adds one unit to x'^. 

The p-value itself is calculated using the x^ distribution with n — 2 degrees of freedom. For semi-central 
collisions a p-value of 0.0009 is found, indicating significant positive It should be noted that 

the most significant data points are at p^^^^ < 60 GeV/c; the results in the range 60 < p^^^^ < 100 
GeV/c are compatible with = 0 (p-value 0.02). In central collisions, a p-value with respect to 
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T 
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Fig. 4: Second-order harmonic coefficient as function a of for 0-5% (a) and 30-50% (b) collision 
centrality. The error bars on the points represent statistical uncertainties, the open and shaded boxes indicate the 
shape and correlated uncertainties (as explained in Sec. 12.51) . 


the hypothesis of v,*’ = 0 of 0.12 is found which indicates that v,” is compatible with 0 within two 


ch jet 


standard deviations. Following the same approach an upper limit of V 2 
same confidence interval. 


ch jet 


0.088 is found within the 


3.1 Comparison to previous measurements and model predictions 

To get a better qualitative understanding of the results, the V 2 of single charged particles 0,121 


and the ATLAS v! 


calo jet 


2 measurement 115811 are shown together with the V 2 


ch jet 


measurement in Fig. [5j The 
ATLAS result is for jets with resolution parameter R = 0.2 within |t] | <2.1 comprising both charged and 
neutral fragments. The event plane angle is measured by the forward calorimeter system at 3.2 < |t]| < 
4.9. Jets are reconstructed by applying the anti-A:T algorithm to calorimeter towers, after which, in an 
iterative procedure, a flow-modulated underlying event energy is subtracted. Each jet is required to lie 
within Arj^ + A(p^ < 0.2 of either a calorimeter cluster of pT > 9 GeV/c or a pj > 10 GeV/c track jet 
with resolution parameter /? = 0.4 built from constituent tracks of pj > 4 GeV/c (the full reconstruction 


procedure can be found in L581 159111. 


It is important to realize that the energy scales of the ATLAS v“*° and ALICE measurements 
are different (as the ALICE jets do not include neutral fragments) which complicates a direct compar¬ 
ison between the two measurements. The central ATLAS results are also reported in 5-10% collision 
centrality. The ALICE and ATLAS measurements are in qualitative agreement, both indicating path- 
length-dependent parton energy loss. Given the uncertainties, the difference in the central values of the 
measurement is not significant. 


Figurel^also shows the V 2 of single charged particles (from i2lL l22ll). which is expected to be mostly 


caused by in-medium energy loss at intermediate and high momenta (px > 5 GeV/c). Even though a 
direct quantitative comparison between and cannot be made as the energy scales for jets and 
single particles are different, the measurements can be compared qualitatively, and it can be seen that for 
central events, the single particle V 2 ”* and V 2 *’ are of similar magnitude and only weakly dependent on 
Px over a large range of px (~ 20 — 50 GeV/c). For non-central collisions (30-50%), the measurements 
of V 2 for single particles and jets are also in qualitative agreement in the px range where the uncertainties 
allow for a comparison. 


Figure [ 6 ] shows the ofR = 0.2 charged jets from the JEWEL Monte Carlo l^l compared to the 


measured JEWEL simulates a parton shower evolution in the presence of a dense QCD medium by 


generating hard scatterings according to a collision geometry from a Glauber II 6 OII density profile. A ID 
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Fig. 5: Elliptic flow coefficient V 2 of charged particles 112 1112211 (red, green) and R = 0.2 full jets (comprising both 
charged and neutral fragments) measured within | rj I <2.1 Q (blue) superimposed on the results from the current 
analysis ofR = 0.2 charged jets v\ measurements, statistical errors are represented by bars and systematic 

uncertainties by shaded or open boxes. Note that the same parton pj corresponds to different single particle, full 
jet and charged jet pj. ATLAS ^‘^'and CMS V 2 from |22, 58| in 30-50 % centrality are the weighted arithmetic 
means of measurements in 10% centrality intervals using the inverse square of statistical uncertainties as weights. 


Bjorken expansion is used to simulate the time evolution of the medium. After radiative and collisional 
energy loss, PYTHIA is used to hadronize the fragments to final state particles. 

The analysis on the JEWEL events is performed with the same jet definition and acceptance criteria 
that are used for the analysis in data, using the symmetry plane 'P 2 from the simulated initial 

geometry as 'Pep, 2- The JEWEL Monte Carlo shows finite significant V 2 *’ in semi-central collisions; in 
central collisions is compatible with zero. The JEWEL result for semi-central 30-50% collisions is 
compatible with the measured values (p-value 0.4 using Eq.[T2with the JEWEL results as hypothesis /r, 
and the quadratic sum of the statistical uncertainties of both datasets as a, in the denominator of the first 
sum of Eq. [T^ . In central JEWEE collisions is consistent with zero, while the measured values are 
compatible with the JEWEE within two standard deviations. It should also be noted that JEWEE 
currently uses an optical Glauber model for the initial state and therefore does not include fluctuations 
in the participant distribution due to the spatial configuration of nuclei in the nucleus. This simplified 


treatment of the overlap geometry may underestimate the 6\\. This comparison of 


ch jet 


in 


JEWEE to experimental data complements earlier studies of the path-length-dependent parton energy 
loss and model predictions for the jet R\\ llH]. 


4 Conclusion 

The azimuthal anisotropy of R = 0.2 charged jet production, quantified as v'^ has been presented 
in central and semi-central collisions. Significant positive is observed in semi-central collisions, 
which indicates that jet suppression is sensitive to the initial geometry of the overlap region of the colli¬ 
sion. This observation can be used to constrain predictions on the path-length dependence of in-medium 
parton energy loss. In central collisions, the central values of the measurement are positive, but the 
uncertainties preclude drawing a strong conclusion on the magnitude of 

The measured V 2 *’ for charged jets is also compared to single particle V 2 from AEICE and CMS and 
^caiojet ATEAS. The measurements cannot be directly compared quantitatively since the energy 
scales are different, but qualitatively, the results agree and indicate a positive V 2 for both charged particles 
and jets to high pj in central and semi-central collisions. This observation indicates that parton energy 
loss is large and that the sensitivity to the collision geometry persists up to high transverse momenta. 
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Fig. 6: vf of R = 0.2 charged jets obtained from the JEWEL Monte Carlo (red) for central (a) and semi-central 

collisions (b) compared to data. JEWEL data points are presented with only statistical uncertainties. 


The JEWEL Monte Carlo predicts sizable for semi-central collisions and very small to zero 
in central events. These predictions are in good agreement with the semi-central measurement. Eor 
central collisions, the JEWEL prediction is below the measurement, but more data would be needed to 
reduce the uncertainties on the measurement sufficiently to constrain the model. 
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